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ABSTRACT 

Spatially resolved multi-wavelength centimeter continuum observations of cool evolved stars can not only constrain the morphology 
of the radio emitting regions, but can also directly probe the mean gas temperature at various depths of the star’s extended atmosphere. 
Here, we use the Very Large Array (VLA) in the A configuration with the Pie Town (PT) Very Long Baseline Array (VLBA) antenna 
to spatially resolve the extended atmosphere of Betelgeuse over multiple epochs at 0.7, 1.3, 2.0, 3.5, and 6.1cm. The extended 
atmosphere deviates from circular symmetry at all wavelengths while at some epochs we find possible evidence for small pockets of 
gas significantly cooler than the mean global temperature. We find no evidence for the recently reported e-MERLIN radio hotspots in 
any of our multi-epoch VLA/PT data, despite having sufficient spatial resolution and sensitivity at short wavelengths, and conclude 
that these radio hotspots are most likely interferometric artefacts. The mean gas temperature of the extended atmosphere ha s a typical 
value of 3000 K at 2 R* and decreases to 1800 K at 6 R*, in broad agreement with the findings of the single epoch study from iLim et alj 
jl998l) . The overall temperature profile of the extended atmosphere between 2R* < r < 6R* can be described by a power law of the 
form TgasCr) « with temporal variability of a few 100 K evident at some epochs. Einally, we present over 12 years of V band 
photometry, part of which overlaps our multi-epoch radio data. We find a correlation between the fractional flux density variability 
at V band with most radio wavelengths. This correlation is likely due to shock waves induced by stellar pulsations, which heat the 
inner atmosphere and ionize the more extended atmosphere through radiative means. Stellar pulsations may play an important role in 
exciting Betelgeuse’s extended atmosphere. 
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1. Introduction 


Our understanding of the mass-loss mechanism in early-M spec¬ 
tral type red supergiants (RSGs) has been inhibited by the lack 
of spatially resolved observations on the scale of the chromo¬ 
sphere and wind acceleration region; the two regions identified 
as the most i mportant for studies of mass -loss mechanisms in 
evolved stars (iHolzer & MacGregoilll985h . Although radiation 
pressure on dust may play a role in driving the wind to its termi¬ 
nal velocity via dust-gas collisions, additional mechanisms are 
required to accelerate the gas out t o distances where s ubstantial 
dust condensation may occur (e.g., iHamer et al.l2009h . Thermal 
free-free centimeter continuum emission arises from the star’s 
extended atmosphere and directly probes the crucial inner spatial 
scales containing the chromosphere and the initial wind acceler¬ 
ation region. Spatially resolved centimeter observations not only 
allow a study of the gas morphology but also enable the mean 
gas temperature to be em pirically derived as a function of depth 
through the atmosphere (iLim et al.lll998ll (hereafter Lh- 98). It 
is therefore a powerful tool to study the details of mass loss in 
RSGs. Despite the value of spatially resolved centimeter contin¬ 
uum observations, very few RSGs have sufficiently large angular 
diameters that they can be resolved with current centimeter in¬ 
terferometers. 


* eamon.ogorman@chalmers.se 


Betelgeuse (a Ori, M2 lab) is t he closest isolated RSG 
(d - 197 + 45 Dc: [jramer et af]i2008h and is therefore the best 
studied of the early-M RSG class. Its has the largest contin¬ 
uum angul ar diameter of all R SGs (0* = 42.49 + 0.06 mas 
at 2.3 pm, lOhnaka et al.l 1201 ih and its extended atmosphere 
has been the focus of many centimeter studies over the past 
few decades. The first detailed study o f Betelgeuse at cen¬ 
timete r wavelengths was carried out by iNewell & Hiellmina 
(Il982h with the Very Large Array (VLA) in its second most 
compact C configuration, at a resolution worse than ~ 1". 
The atmosphere was unresolved but the radio emission was 
interpreted as chromospheric in origin and extending from 
1 to 4R*. This wa s in broad agreement w ith subsequent 
Alfven wave models (IHartmann & Avretliri984h and later Hub¬ 
ble Spa ce Telescope (HST) spatially resolved ultraviolet obser- 
vations (iGilliland & Dupr^ll996l:rUitenbroek et al.lll998l) . Spa¬ 
tially resolved VLA plus Multi-Element Radio Linked Interfer¬ 
ometer Network (MERLIN) observations at 6 cm a lso confirmed 
the ex tended nature of the radio emitting region (ISkinner et al.l 
l1997h . 

Lh- 98 used the VLA in its most extended A configuration 
to resolve Betelgeuse’s atmosphere at 0.7 cm (with 40 mas res¬ 
olution) and partially resolve it at 1.3, 2.0, 3.5, and 6 cm (with 
resolution between 90 and 330 mas). Because the radio emis¬ 
sion is thermal and optically thick, they were able to derive the 
mean gas temperature as a function of radius, and discovered 
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that the temperature never reached chromospheric values (i.e., 
~ 4000 - 8000 K) but decreased steadily from ~3450K at 2 R*, 
to ~1370K at 7R*. They also detected an asymmetry in their 
0.7 cm image which they attributed to the action of a large con¬ 
vection cell. To reconcile their results with the extended ultravio¬ 
let observations, which probed the same regions, they concluded 
that the inner atmosphere must be inhomogeneous to accom¬ 
modate the hot chromospheric plasma, but th at the cooler gas 
must b e 3 orders of magnitude more abundant. iHamer & BrownI 
(I2OO6I1 used observations of the chromospheric electron density 
tracer C II] 42325 A to confirm this low filling factor for the 
chromospheric gas. _ 

Recently, an unexpected discovery by [Richards et al.l (1201 3h 
(hereafter Rh- 13) with e-MERLIN at 80 mas x 60 mas resolution, 
revealed that two unresolved radio features contribute about 75% 
of the total flux density at 5.2 cm. The two “radio hotspots” are 
separated by 90 + 10 mas (i.e., ~ 4 R*) and have inferred bright¬ 
ness temperatures Tt, = 3800 + 500 K and Tj, - 5400 + 600 K, 
with the later value signifi cantly above the pho tospheric temper¬ 
ature (Teff = 3690 + 54 K, lOhnaka et ai1l201 Ih . Indeed, because 
the features are unresolved, and assuming the emission is ther¬ 
mal and optically thick, then these inferred brightness tempera¬ 
tures will be lower limits for the gas temperature. The potential 
importance of these radio hotspots for understanding the wind 
acceleration region of RSGs has been the motivation for this 
paper. In the following sections we present multi-epoch multi¬ 
wavelength radio observations of Betelgeuse, taken ~ 10 years 
prior to these e-MERLIN observations. At short wavelengths, 
our data have comparable or superior spatial resolution to the 
e-MERLIN data and so the goal of this paper is to search for 
signatures of hotspot features to improve our understanding of 
Betelgeuse’s wind acceleration region. 

2. Observations and data reduction 

Betelgeuse was observed with the NRACQ VLA in the A con¬ 
figuration with the Pie Town (PT) Very Long Baseline Array 
(VLBA) antenna at multiple epochs between December 2000 
and October 2004. The PT antenna is located approximately 
52 km from the center of the VLA, and its inclusion to the A 
co nfiguration enhances the east-west resolution by a factor of 
2 (lUlvestad et alJlT99^ . The dates and corresponding combina¬ 
tions of wavelengths used in our observations are summarized 
in the first two columns of Table [T] Our two multi-wavelength 
data sets in 2002 (Program; AH0778) were taken only ~ 2 
months apart and consist of all wavelength bands between 1.3 
and 20.5 cm. We observed the star again under program AH0824 
in 2003 and 2004 at all wavelength bands between 0.7 and 
20.5 cm. We also obtained VLA A configuration plus PT antenna 
data from the NRAO data archive (Program: AL0525) at 0.7 and 
1.3 cm which were taken in late 2000 and early 2001, respec¬ 
tively. 

We used the fast switching phase calibration technique at 
0.7, 1.3 , and 2.0 cm to com pensate for tropospheric phase vari¬ 
ations dCarilli et al.lll99dl . switching between observations of 
the point source phase calibrator 0532H-075 (located 5.6° from 
Betelgeuse) and Betelgeuse, with cycle times varying between 2 
and 3 min. The archival data from 2000/2001 also implemented 
the fast switching technique at these wavelengths but used the 
point source phase calibrator 0552H-032, which was ~ 10 times 

' The National Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Foundation operated under cooperative agreement by As¬ 
sociated Universities, Inc. 


fainter at 0.7 cm but 4.2° from Betelgeuse. We also interleaved 
observations of 0532-1-075 to calibrate phase variations at longer 
wavelengths, but implemented longer cycle times of 10 - 20 min. 
Absolute flux density calibration was obtained from 3C48 at 
6.1 cm, and 3C286 at all other wavelengths. The absolute flux 
density uncertainty at all epochs is estimated to range from 3% 
to 10% between the longest to shortest wavelengths, respec¬ 
tively. The flux density uncertainties quoted in the remainder 
of this paper do not include this systematic error but represent 
the +lcrrms random error. The raw visibilities were flagged and 
calibrated according to standard VLA continuum data reduction 
proced ures using the Astronomical Image Processing System 
(AIPS: lGreisenlll990h . The calibrated visibilities were imaged in 
both AI PS and the Common A stronomical Software Application 
(CASA: liVIcMullin et al.ll2007ll package. 

2.1. Imaging 

The calibrated visibilities at all wavelengths were imaged us¬ 
ing CASA’s clean task while applying uniform weighting. Ad¬ 
ditional images were also created a t 0.7 and 1.3 c m which im¬ 
plemented both Briggs weighting (iBriggsl ri995h . and natural 
weighting with a restoring beam size equal to that obtained us¬ 
ing uniform weighting, as used in Lh- 98. These additional im¬ 
ages were used to check if any asymmetries produced by one 
weighting scheme could be reproduced by another. 

2.2. Calibrated visibilities 

Betelgeuse was observed offset by a few restoring beams from 
the phase reference center to avoid any possible spurious arte¬ 
facts which might accumulate there. One such artefact was no¬ 
ticed in the 2.0 cm image from April 2002 but our target was 
sufficiently far away and was not contaminated. We fitted uni¬ 
form elliptical Gaussian brightness distributions to the source in 
each image and measured the position of the peak emission. We 
then used CAS As fixvis task to shift the visibilities such that this 
peak emission position was placed at the phase reference center. 
This shift enabled us to extract the azimuthally averaged com¬ 
plex visibilities discussed later in Sectio n l3.2l The Python-based 
task uvmultifit (iMartf- Vidal et al.ll2014ll was then used to fit uni¬ 
form intensity disks to the calibrated visibilities, with the major 
axis (0maj), axis ratio (^min/fi’maj, where ^min is the minor axis), 
position angle (PA), and flux density {Fy) left as free parame¬ 
ters. 

3. Results 

3.1. Radio images 

Radio images were produced at all wavelengths using uniform 
weighting to utilize the full capabilities of the PT baselines and 
achieve maximum spatial resolution. The corresponding restor¬ 
ing beam sizes and image root mean square (rms) noise levels are 
provided in columns 3 and 4 of Table [T] The images were gener¬ 
ally of high quality with the source clearly resolved at the shorter 
wavelengths. There were no obvious large scale asymmetries 
similar to those detected by R-i-13 present in any of the images at 
any wavelength. To investigate the presence of possible smaller 
scale asymmetries, analysis was carried out on the visibilities 
as the source to beam size ratio was small and this is discussed 
in Section 13.21 The 0.7 cm images were the only ones to con¬ 
tain any obvious small scale (i.e., < 40 mas) asymmetries. How¬ 
ever, when these data sets were imaged using the different im- 
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Table 1. Multi-epoch VLA A configuration plus PT antenna observations of Betelgeuse. 


Date 

A 

Restoring beam 

^rms 

^maj 

^min 

PA 

Fv 

Tb 


(cm) 

(mas X mas, °) 

(mJybeam^') 

(mas) 

(mas) 

(°) 

(mJy) 

(K) 

2004 Oct 21,30 

0.7 

39x26,42 

0.37 

99 + 3 

92 + 3 

92 + 20 

28.68 + 0.53 

2940 + 140 


1.3 

80 X 42,32 

0.09 

127 + 2 

114 + 2 

85 + 5 

13.83 + 0.10 

3160+80 


2.0 

121 x91,54 

0.08 

158 + 6 

152+ 10 

132 + 59 

7.23 + 0.15 

2360+ 190 


3.5 

208 X 126,41 

0.02 

218 + 6 

182 + 9 

162 + 7 

3.34 + 0.03 

2140+120 


6.1 

377x264,46 

0.02 

315 + 29 

187 + 40 

173 + 9 

1.55 + 0.04 

1920 + 450 


20.5 

1262x889,45 

0.03 

< 889 

< 889 


0.25 + 0.03 

>260 

2003 Aug 10,12 

0.7 

40x27,47 

0.46 

103+4 

92 + 5 

104+16 

28.05 + 0.84 

2760 + 200 


1.3 

80 X 42,32 

0.17 

122 + 5 

114 + 6 

8 + 24 

11.20 + 0.24 

2670+ 190 


2.0 

119x96,62 

0.10 

132+ 10 

116 + 9 

11+27 

5.88 + 0.17 

3010 + 340 


3.5 

204 X 139,46 

0.03 

193 + 7 

140+11 

152 + 7 

2.80 + 0.04 

2630 + 230 


6.1 

378x297,69 

0.03 

247 + 47 

179 + 61 

169 + 28 

1.22 + 0.04 

2010 + 780 


20.5 

1247x931,49 

0.04 

< 931 

<931 


0.26 + 0.03 

>250 

2002 Apr 12,13 

1.3 

91 X 59,45 

0.18 

134 + 9 

102 + 6 

36+10 

8.96 + 0.24 

2170 + 200 


2.0 

131 X 98,60 

0.39 

166+16 

104+ 13 

41 + 11 

5.32 + 0.24 

2420 + 400 


3.5 

224 X 155,65 

0.03 

234 + 9 

171 + 10 

41+7 

2.66 + 0.04 

1690+ 120 


20.5 

1398 X 1146,55 

0.06 

< 1146 

< 1146 


0.38 + 0.06 

>240 

2002 Feb 17,18 

1.3 

83 X 48,30 

0.14 

120 + 4 

109 + 3 

31 + 13 

10.87 + 0.17 

2750+ 130 


2.0 

128x90,51 

0.11 

148+11 

125 + 18 

142 + 26 

5.38 + 0.22 

2280 + 380 


3.5 

200 X 135,43 

0.03 

222+8 

166+13 

155 + 7 

2.85 + 0.04 

1960+ 170 


20.5 

1312x951,55 

0.05 

< 951 

<951 


0.30 + 0.05 

>270 

2001 Jan 02 

1.3 

78 X 42,35 

0.08 

124 + 2 

115 + 1 

40 + 8 

12.58 + 0.08 

2920 + 60 

2000 Dec 23 

0.7 

44 X 20,29 

0.18 

98 + 2 

91 + 1 

178 + 8 

29.02 + 0.30 

3040 + 80 


Notes. The restoring beam properties (major axis, minor axis, PA) and image rms noise values are taken from the uniformly weighted radio images 
that include the PT antenna baselines. The major and minor axis of the resolved stellar radio disk and the position angle PA (measured 
east of north), and the total flux density Fy, are all derived from the best-fit uniform intensity disk models using uvmultifit. The atmosphere at 
20.5 cm was unresolved so we define the minor axis of the restoring beam as the upper limit to its diameter. The brightness temperatures, T*, have 
been derived using Equation[T]in Section [T4l 


age weighting schemes specified in Section lZTl the morphology 
was inconsistent. For example, the images created using standard 
uniform weighing were morphologically different to the images 
created using natural weighting with a restoring beam of size 
corresponding to that obtained with uniform weighting. These 
inconsistencies were probably due to poor phase stability result¬ 
ing from rapid tropospheric phase fluctuations which our cycle 
times were not short enough to compensate for. To investigate 
this further, we removed the PT antenna baselines and again im¬ 
aged the data using the two aforementioned imaging techniques. 
We found that two 0.7 cm data sets (i.e., from 2000 and 2004) 
produced images that had consistent morphology when any of 
the imaging techniques were applied and showed no asymme¬ 
tries. The observing logs of the other data sets report poorer 
weather conditions in agreement with the conclusions from our 
imaging analysis. 

Our reliable 0.7 cm images are plotted in Figure [T] alongside 
the 0.7 cm image from Lh- 98 for comparison. In these data sets, 
the cycle time was adequate to obtain reasonable phase trans¬ 
fer from the calibrator to the source for the VLA baselines but 
shorter cycle times were probably required for the other data sets 
due to poorer weather. We also independently reduced the L-i-98 
0.7 cm data set and were able to reproduce the notable asym¬ 
metry to the East, irrespective of the applied image weighting. 
Unlike the L-i-98 image, no large scale asymmetries are present 
in either of our two 0.7 cm images as shown in Figure [1] and so 
the asymmetry reported by L-i-98 was not a large scale (~ 2 R*) 
stable feature. 


3.2. Radio diameters 

It can be seen in Figure [1] that the size of the source at 0.7 cm 
is only a few beam diameters across. The source to beam ratio 
gets progressively smaller at longer wavelengths which makes 
fitting models to the images unreliable. In this case it is best 
to obtain estimates of the source size and flux density by di¬ 
rectly fitting to the calibrated complex visibilities. In Figure |2] 
we plot the real component of the complex visibility against the 
projected baseline length in wavelengths, B^, for the 2004 epoch 
at all wavelengths between 0.7 - 6.1 cm. The data has been az- 
imuthally averaged into uniformly sized bins over large portions 
of the unevenly sampled u - v plane, where u and v are the spa¬ 
tial coordinates of the projected baseline in wavelengths, and 
+ v^. It is clear that the source is resolved at all wave¬ 
lengths shown, i.e., we have obtained visibilities beyond to the 
first null. The VLA observations of Lh- 98 only fully resolved 
Betelgeuse’s atmosphere at 0.7 cm at one epoch, but here we 
fully resolve the atmosphere at all wavelengths between 0.7 and 
6.1 cm over multiple epochs thanks to the inclusion of the PT an¬ 
tenna. The black data points in Figure|2]represent the real visibil¬ 
ity measurements from the VLA antenna baselines only, the red 
data points include both VLA and PT antenna baselines, while 
the green data points include only PT antenna baselines. This 
highlights the importance of the inclusion of the PT antenna, es¬ 
pecially at the longer wavelengths where it provides many mea¬ 
surements close to the first null, which helps to constrain the size 
of the atmosphere. The corresponding binned imaginary compo¬ 
nent of the complex visibilities are not shown in Figure |2] but 
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Fig. 1. VLA A configuration images of Betelgeuse at 0.7 cm at three different epochs. These images have been created by naturally weighting 
the visibilities and applying a restoring beam corresponding to the size obtained from uniform weighting. The restoring beam is located in the 
bottom left corner of each panel while the red circle is the approximate location and size of the optical photosphere. The contour levels are plotted 
at 10%, 20%,..., 90%, and 99% of the peak flux density. From left to right the beam sizes are 40" x 40", 45" x 41", and 48" x 40". Left panel: 
The L+98 image shows an asymmetry to the east of the pho tosphere which is here assumed to be coincident with the intensity-weighted centre of 
the radio disk. A photospheric angular diameter of 45 mas iPvck et alJll99fih was assumed by the authors. Middle and right panels: Our images 
show no large scale asymmetries. The posit ion of the photosph ere is here assumed to be located at the peak of the emission and we adopt the 
photospheric angular diameter of 42.49 mas dOhnaka et alj20lll) 


usually had zero values at all baselines, indicative of no large 
scale asymmetric emission. 

Uniform intensity disk models were fitted to the complex vis¬ 
ibilities as outlined in Section The major and minor axes, 
0niaj and flmin, and PAs of the best fits are listed in columns 5, 
6, and 7 of Table [1] In Figure |2] we include these theoretical 
visibility models, assuming of a uniform intensity circular disk 
with an effective angular radius defined as = 0.5 
These models provide reasonable though not perfect fits to the 
azimuthally averaged real visibilities at all epochs. Beyond the 
first null we can see that there is an even greater departure from 
the uniform intensity disk model at short wavelengths and this 
departure is evident at other epochs too. Limb darkening and 
brightening effects, and a weak unresolved feature located away 
from the phase center could potentially manifest themselves in 
the second lobe. Limb brightening could possibly be caused by 
the temperature rise from the photosphere to the 1.4 t emper - 
ature peak in the semi-empirical model of iHamer et akl (1200 Ih . 
We find no clear signature in the imaginary component of the 
complex visibilities or phases at the corresponding baselines for 
either limb darkening or an additional unresolved weak feature. 
Indeed, we find that the data provides better fits to the uniform 
intensity disk model beyond the first null for epochs which have 
the best u - v sampling and so the consistently poor fit beyond 
the first null in Figure |2] may well be due to the limited u - v 
sampling. 

To further investigate the possible presence of small scale 
asymmetries in the radio emission morphology, we subtracted 
the best fit uniform disk model with parameters listed in Table[T] 
from the corresponding visibilities for each data set. After imag¬ 
ing, most of these uniform disk subtracted data sets consisted of 
only noise, but a few contained one or two weak but significant 
(i.e., 3-6 cTjjns) negative flux density point sources. To derive 
reliable properties for these point sources and to rule out image 
deconvolution artefacts, we went back and simultaneously fitted 
a uniform disk and a delta function to the original calibrated vis¬ 


ibilities. For five datasets we found that a uniform disk and a 
negative flux density delta function offset from disk center pro¬ 
duced a slightly better fit than just a single uniform disk. The 
parameters from these fits are listed in Table lA.ll The bright¬ 
ness temperatures of the point sources were found to vary be¬ 
tween ~ 300 - 800 K lower than the simultaneously fit uniform 
disk brightness temperature. Their unresolved nature means that 
these brightness temperatures are upper limits for the gas tem¬ 
perature, for thermal and optically thick emission. We interpret 
these negative point sources as pockets in the extended atmo¬ 
sphere where the local gas temperature is lower than the mean 
gas temperature and highlight the non-uniform nature of Betel- 
geuse’s extended atmosphere. 

The centim eter continuum free-free opacity varies as 

oc ' (e.g.. iRvbicki & Lightmanlll979l) and so our multi¬ 
wavelength radio observations probe different layers of Betel- 
geuse’s stellar atmosphere along the line of sight, with the longer 
wavele ngths probing emis sion further away from the photo¬ 
sphere. iH^^ereF^dTOOil) showed that the dominant source of 
this opacity is from the free-free interactions of free electrons 
from photoionized metals with neutral hydrogen. We find the 
source size gets progressively larger with increasing wavelength 
which implies that the radio surfaces are optically thick, and is in 
agreement with the findings ofL-i-98. This is clearly illustrated in 
Figure|3]where the values of the effective angular radius are plot¬ 
ted against wavelength over all epochs between 2000 and 2004. 
We find that the power law (peg - (54.3 + ® mas can 

describe the size of the radio emitting region well at all wave¬ 
lengths between 0.7 and 6.1 cm as shown by the dashed line in 
Figure [3] It is clear from Figure [3] that millimeter and submil¬ 
limeter spatially resolved observations will be required to study 
the morphology of the inner atmosphere within 2 R*. The Ata¬ 
cama Large Millimeter/submillimeter Array (ALMA) will have 
such capabilities in its more extended array configurations at all 
wavelengths between Band 2 and Band 10, as shown in Figure 

13 
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B,(X) 


Fig. 2. Real component of the complex visihility plotted against the 
projected baseline length in wavelengths, Bx, for the 2004 epoch at all 
wavelengths between 0.7 - 6.1 cm. The black data points contain only 
VLA antenna baselines, the red data points contain both VLA and PT 
antenna baselines, while the green data points contain only PT antenna 
baselines. The dashed blue line represents the theoretical visibilities of a 
uniform intensity circular disk with an ejfective angular radius defined 
as (^eff “ 0.5 '^Umaj ^min• 


We can also investigate if there is any variability in the size 
of the extended atmosphere over time at the various centimeter 
wavelengths. At 0.7 and 1.3 cm there is no significant change 
in 0efF at all epochs with A^eft ^ 3%, where A^eff is the max¬ 
imum deviation about the mean. There is more variability in 
0eff at longer wavelengths which sample the more extended at¬ 
mosphere, with Aipsff ~ 15% at 2.0 and 3.5 cm. We only have 
two measurements at 6.1 cm but the associated error bars are too 
large to detect any significant variability in 0eff- 

3.3. Flux density variability 

Betelgeuse’s atmosphere is optically thick at centimeter wave¬ 
lengths and so the thermal continuum flux density can be ap¬ 
proximated as emanating from an opaque disk with Fy oc 
Tgasiplfi- The relationship tells us that variations in the radio 
flux density are caused by either changes in the mean gas tem¬ 
perature or changes in the size of the extended atmosphere. 


7 
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3 „ 

2 i 
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0.1 1.0 10.0 
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Fig. 3. The effective angular radius of Betelgeuse’s extended atmo¬ 
sphere as a function of wavelength. It can be described by the power 
law, 0e(f = 54.3/1^^ cm, between 0.7 and 6.1cm. Multi-wavelength 
extended ALMA configuration observations will be required to study 
the morphology of t he radio emission w ithin 2 R*. The molecular shell 
(MO Lsphere) from iPerrin et all ( l2007h and the photospheric radius 
from lOhnaka et all (1201 ih are also shown. 
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iNewell & HiellminS (Il982l) found that multi-wavelength radio 
flux density measurements from the period 1972- 1981 could be 
descri bed by the radio sp ectrum Fy - (0.23 + 0.01)Vg^^^'^° °'mJy, 
while iDrake et^ (119921) reported an almost identical spectrum 
for data spanning 1986 - 1990. In Figure |4] we plot these ra¬ 
dio spectra along with our VLA/PT flux density values and the 
values from L-i-98. It can be seen that all the data from 1996 
onwar ds lie below the s pectru m of both iNewell & Hiellming 
(1^9^) and iDrake et al.l il992h as initially reported in Hamer 
( 2013 ) and are better fitted by the spectrum Fy — (0.18 + 
0.01)vQ^^*°*^^mJy. Therefore, even though the spectral index, a 
(where Fy oc v“), is unchanged, the flux density appears to have 
reduced by ~ 20% over the last few decades. 

The previous flux density measurements prior to this work 
were obtained with sin gle dish instrument s and with the VLA in 
various configurations. IPrake et al.l (Il992h found that their 2 cm 
flux densities were systematically lower in the more extended 
VLA configurations than in the more compact ones, but this was 
probably a combined result of using long cycle times (up to 17 
minutes in duration) and taking the peak of their Gaussian fits 
for the flux density of a resolved or partially resolved source. 
Our shorter cycle times (e.g., 3 minutes at 2 cm) along with our 
method of fitting directly to the visibilities to extract the flux 
density will avoid these problems and so different VLA config¬ 
urations cannot account for the recent ~ 20% reduction in flux 
density. It is still conceivable that at the very high frequencies, 
our cycle times might still have been less than the coherence 
time, and could have resulted in a loss of correlated amplitude. 
However, at longer wavelengths (i.e., > 2.0 cm) we still And the 
flux density to be consistently lower than historical values and so 
poor coherence cannot account for our lower flux density mea¬ 
surements. Finally, systematic errors are unlikely to account for 
the 20% reduction in flux density. The absolute flux density un¬ 
certainty is only 3 - 5% at the longer wavelengths, while the 
main flux dens ity calibrator 3C286, v aries by only a few percent 
every decade (IPerlev & Butleiir2013h . This would suggest that 
there has been either a global reduction in the mean gas temper¬ 
ature in Betelgeuse’s extended atmosphere or a global reduction 
in the size of the radio emitting region in the last few decades. 
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Fig. 4. Radio spectra for Betelgeuse at VLA wavelengths spanning 
1972 — 2012. The continuo us line represents the radio spectrum from 
iNewell & Hiellmiri3 lll982h and includes measurements u p to 90GHz. 
The d ashed line is the almost identical spectrum from [Drake et'H] 
(Il992h . The dash-dotted line is the radio spectrum that best fits the more 
recent data between 1996-2004 and is described by Fy = 0.18 mjy v* . 
The flux density from the e-ME RLIN 5.2 cm observ ation and a 100 GHz 
flux density measurement from iHarper et alj l l2009h are also shown. 


During our multiple epochs of observations, we detect low 
level flux density variability (< +25%) at all wavelengths ex¬ 
cept at 0.7 cm, which shows no significant variability. This is 
in excellent agreement with the findings of [Drake et al.l (Il992h 
who monitored Betelgeuse regularly over a four year period at 
2.0, 3.5, and 6.1 c m and reported similar levels of variability. 
[Drake et"^ ([1992h also found no correlation between variabil¬ 
ity at the three wavelengths sampled and noted significant vari¬ 
ability on time scales of < 1 month. Our multi-wavelength data 
is sparsely sampled over the ~ 4 years but we do And signifi¬ 
cant variability at 1.3 cm in just two months. We also And that 
the flux density increases between our 2003 and 2004 epochs 
at 2.0, 3.5, and 6.1cm and this increase coincides with an in¬ 
crease in the corresponding effective radii. However, a change in 
flux density does not always coincide with a change in the at¬ 
mospheres size. For example, the 1.3 cm flux density variability 
between the two sampled months of 2002 occ urs in spite of the 
fact th at the (p^g remains essentially constant. [Bookbinder et al.[ 
([1987h used time scale arguments to suggest that this short term 
variability in the radio flux density could not be caused by large 
scale processes such as changes in the mass-loss rate or global 
pulsa tions because these dyna mical time scales are much too 
long. [Bookbinder et al.[ ([l987l) suggested this short term vari- 
ability could be due to the molecular catastrophe scenario of 
[Muchmore et al.[ ([1987h which predicts that global depletions of 
ions into molecules could rapidly chan ge the radio flux density . 
This mechanism has been explored by Hamer & Linskvl ([200 ih 
who found that either incomplete or large scale formation of CO, 
SiO, and H 2 could be responsible for the observed radio vari¬ 
ability triggered by rapid changes in the far ultraviolet radiation 
held. We propose an alternative mechanism for the flux density 
variability in Section [4^ 


3.4. Temperature profile of the extended atmosphere 

Spatially resolving the extended atmosphere at multi¬ 
wavelengths also allows the brightness temperature to be 
calculated as a function of depth along the line of sight. If we 
assume the radio emission is emanating from a uniform intensity 
disk with major and minor axis, ^maj and respectively, then 
the brightness temperature can be written as 


Tb = 1.96 X 10^ 




1 ^maj ^min \ 

ImJyj 

Icm/ ' 

[ mas^ ) 


-1 


K. 


( 1 ) 


Since the radio emission from Betelgeuse is thermal and op¬ 
tically thick, the brightness temperature is a measure of the 
mean gas temperature, Tgas, and so multi-wavelength spatially 
resolved radio observations allow the radial mean gas tempera¬ 
ture profile to be constructed. 

In the Anal column of Table [1] we list the brightness temper¬ 
ature at each wavelength and epoch which have been derived 
using Equation [1] These values are plotted against the effec¬ 
tive angular radius, or distance from the photosphere, in Fig¬ 
ure |5] alo ng with the values from L+98 and the semi-empirical 
model of lHamer et al.[ ([200 Ih . The gas temperature profiles con¬ 
structed for all our epochs between 2000 and 2004 are in good 
agreement with both the values derived from the si ngle epoch of 
L+98 and the subsequent semi-empirical model of [Harper et al.[ 
([200 Ih . Our data implies that the mean gas temperature has a 
typical value of 3000 K at 27?*, and decreases to 1800K at 67?*. 
We do not resolve the atmosphere at 20.5 cm and so our data 
does not probe the gas temperature beyond 67?*. We also note 
that shorter wavelength (i.e., < 0.7 cm) spatially resolved obser¬ 
vations will be required to probe the mean gas temperature of the 
atmosphere inside 2 R*. It is within thi s region that cool molec¬ 
ular layers of CO, SiO, an d H 2 O exist ([Tsuii[[20()0t [Perrin et"^ 
\20ffk [Qhnaka et al.[[20li1) and nascent molecule formation may 
take place, as shown in Figured 

We also detect temporal variability in the gas temperature at 
certain regions of the atmosphere. At 4.7 7?* we And the tem¬ 
perature to change by 400 K over 2.5 yr while closer in towards 
the star at 2.87?*, we And the temperature to change by almost 
600 K in under two months. These changes are significant, even 
when we allow for the absolute flux density uncertainties. Never¬ 
theless, most values of the gas temperature are consistent within 
the random uncertainties implying that the mean gas temperature 
profile is generally stable with occasional fluctuations on the or¬ 
der of a few hundred Kelvin. Fitting a power-law to all the values 
displayed in Figure |5] yields an empirically derived temperature 
profile of the form TgasCr) oc in good agreement with 

the power-law f it of T^^d.r') or r_^^o the semi-empirical temper¬ 
ature profile of [Harper et akl ([200 1[) . This temperature profile is 
less steep than the r ' falloff derived for the red giant Arcturus 
(a Boo: K2 III) based on mo del dependent multi-wa velength un¬ 
resolved centimeter studies ([Q’Gorman et al.|[2013h . 


4. Discussion 

4.1. Where are the e-MERLIN radio hotspots? 

The multi-epoch spatially resolved VLA/PT data presented in 
this paper contain no large scale asymmetries in the radio emis¬ 
sion morphology, let alone radio hotspots, at any wavelength. 
In Figure |6] we plot the e-MERLIN 5.2 cm image adapted from 
R+13 which clearly shows the two radio hotspots, separated by 
90 mas and having a PA of 110° along their line of centers. For 
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‘t>eff (mas) 

Fig. 5. The mean gas temperature profile of Betelgeuse’s extended at¬ 
mosphere revealed hy multi-wavelength spatially resolved VLA and 
VLA/PT observations spanning 1996 - 2004. The red filled circles are 
from the work of L-l-98 while the dotted line is the predicted gas tem per- 
ature profile from the semi-empirical model of lHamer et al.l (1200 ih . All 
other data points are from this work. The error bars represent the ± liTnns 
level, except the 0.7 cm error bars of L-(-98 which include a ±20 % 
uncertainty in the absolute flux density. The vertical and horizontal 
dashed lines represent the photospheric radius of 21.24 ma s and the pho- 
tosph eric effective temperature of 3690 K, respectively dOhnaka et al.l 

l20lTh 


comparison, we also include the restoring beams used to cre¬ 
ate the radio images from our 2004 epoch at 0.7 cm, 1.3 cm, 
and 2.0 cm. The absence of these radio hotspots in our VLA/PT 
data at 20.5 cm, 6.1 cm, and probably 3.5 cm could be explained 
by insufficient spatial resolution. However, our multi-epoch ob¬ 
servations at both 0.7 cm and 1.3 cm have more than adequate 
spatial resolution to resolve these features. Moreover, the spa¬ 
tial resolution of our 2.0 cm VLA/PT images along the direction 
of the hotspots’s line of centres is ~ 95 mas and should also be 
sufficient to resolve some signatures of the reported hotspots. 
Since we had the capabilities to spatially resolve the reported ra¬ 
dio hotspots at three wavelengths over multiple epochs, we now 
examine why they have no signature in our VLA/PT data. 

The unresolved radio hotspots detected with e-MERLIN 
at 5.2 cm (i.e., C band) had flux density values of 0.53 and 
0.79 mJybeam”* with a systematic flux density uncertainty of 
10%. In Section [33] we demonstrated that a spectral index of 
a - 1.33 is appropriate for the total flux density of Betelgeuse at 
centimeter wavelengths over the three decades prior to 2004. If 
we assume that the radio hotspots obey this spectral index then 
we should have detected even the weakest radio hotspot at the 
26crrnis level at 1.3 cm and at the level at 0.7 cm, in our 

uniformly weighted images. Assuming the emission to be ther¬ 
mal and taking the extreme case whereby the radio hotspots are 
optically thin where a — -0.1, we would only have just been ca¬ 
pable of detecting the weakest radio hotspot at the 3o-^^ level at 
1.3 cm while it would not have been detectable at 0.7 cm. How¬ 
ever, it would be highly unlikely for the radio hotspots to be 
optically thin, because in such a case Tgas >> T/,, and the gas 
would be fully ionized and thus many orders of magnitude more 
opaque. A final scenario is the possibility that the emission from 
one or both radio hotspots is non-thermal and so the spectral 
index could be very negative and would explain their absence in 



Fig. 6. Contours represent the e-MERLIN 5.2 cm image of Betelgeuse 
obtained in 2012 with levels set at (5, 10, 15, ....,70) x <Trn,s (where 
CTnns = 9yuJy). The e-MERLIN restoring beam is shown in the bottom 
left corner, while the restoring beams of our 2004 epoch data at 0.7 cm, 
1.3 cm, and 2.0 cm are shown at the top left for comparison. The red 
filled circle is the scaled p hotosphere at the 20 12 position predicted by 
the astrometric solution of lHamer et all(l2008h . The blue, green, and or¬ 
ange ellipses are the derived sizes of the atmosphere at 0.7 cm, 1.3 cm, 
and 2.0 cm, respectively, for the 2004 epoch. The combined positional 
uncertainty is 28 mas. 


our high resolution data. However, our 20.5 cm flux density mea¬ 
surements show no major rise above the expected spectral index 
value of 1.33 and so we can rule out any significant non-thermal 
emission. The absence of the radio hotspots in our multi-epoch 
high spatial resolution data sets cannot be explained by insuffi¬ 
cient sensitivity. 

HST observations provide us with additional spatial in¬ 
formation about Betelgeuse’s extended atmosphere. Ultravio¬ 
let Faint Object Camera images show a non-uniform brightness 
distribution and have appr oximately twice the optical diameter 
(iGilliland & DuDre3ll996li . They therefore probe the same re¬ 
gion as our 0.7 cm radio data presented above. However, no re¬ 
ports from HST imaging have reveal ed features separated by 
90mas (e.g.. iDunree & Stefanikll^l3h . While collisionally ex¬ 
cited optically thin C II] 2325A emission is observed out to 
+ ~ 4R*, no strong ex cess emission was detec ted far from the 
ultraviolet disk centre (iHamer & Brownll2006h . HST spatially 
scanned spectra reveal radial velocity shifts in ultraviolet spec¬ 
tral features for certain scan angles and not in others. If we inter¬ 
pret the absence of radial velocity shifts to be the scans along the 
stellar rotation axis, then the PA of the rotation a l axis is probably 
~ 65° (see the discussion in iHamer & Browril (l2006l) and refer¬ 
ences therein). This PA does not align with the 110° PA of the 
e-MERLIN radio hotspots line of center and so they can not be 
associated with the rotational pole of the star, where hot spots 
might be expected du e to rotation induced asymmetries (e.g., 
lUitenbroek et al.l[T998l) . 

The data presented in this paper were obtained between 8 and 
12 years prior to when the e-MERLIN radio hotspots were re¬ 
ported. We have also re-reduced the multi-wavelength VLA data 
from L+98 and can confirm that no radio hotspots were present. 
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The possibility that the extended atmosphere could have under¬ 
gone such an enormous morphological and thermal reconfigura¬ 
tion in less than a decade after our VLA/PT observations seems 
remote. Such a change would be in stark contrast to the ‘radio 
quiet period’ of both our VLA/PT observations and the VLA 
observations of L-i-98 which together span a total of ~ 8 years 
and show the atmosphere to be cool, continuously more extended 
at longer wavelengths, with little or no large scale asymmetries 
in the morphology. If one (or both) of the radio hotspots were 
indeed new features, then one would presumably expect a sig¬ 
nificant increase in the total flux density, but such an increase 
is in disagreement with the historical flux density measurements 
described in Section [33] 

With such a large amount of multi-epoch multi-wavelength 
data at our disposal showing no trace of the reported e-MERLIN 
radio hotspots, one of us (AMSR) has now re-reduced the e- 
MERLIN data. It has recently been discovered that the position 
of the Cambridge telescope receiver axis was not entered cor¬ 
rectly in the e-MERLIN data during the 2012 Betelgeuse obser¬ 
vations (lBeswickll2015h . This receiver axis offset has now been 
corrected off-line and preliminary re-imaging at high spatial res¬ 
olution shows no trace of the previously reported radio hotspots. 
An examination of the calibration procedure is currently under 
way to assess the impact of the receiver axis offset correction on 
the final e-MERLIN data set. Irrespective of this, we can con¬ 
clude that the radio hotspots as reported in R-h 13 are most likely 
interferometric artefacts. 

4.2. Radio flux density versus optical photometry 

We have also obtained optical V band photometry of Betelgeuse 
(V = 0.42, B-V = 1.85) as part of our temporal evolution study 
of its extended atmosphere. This allows us to search for rela¬ 
tionships between changes close to the stellar photosphere and 
changes in the extended atmosphere. The V band observations 
were taken with 20 and 28 cm Schmidt-Cassegrain telescopes 
near Villanova University. An uncooled Optec, Inc. SSP-3 solid 
state photometer with a silicon red-sensitive PIN-photodiode de¬ 
tector was used to record raw photometer counts through a V 
band filter centered on 550 nm. The filter had a similar spectral 
response to the standard Johnson V band filter. The photome¬ 
try itself was conducted differentially with respect to the KO Illb 
comparison star HD 37160 (V = 4.09, B-V = 0.95). The com¬ 
parison star was further observed differentially with respect to 
the B2III check star HD 35468 (V = 1.64, B-V = -0.22) to 
check the non-variability of the comparison star. Apart from ran¬ 
dom magnitude scatter, there was no variation between the com¬ 
parison and check stars. 

In the bottom panel of Pigure|7]we plot our V band photom¬ 
etry of Betelgeuse spanning from September 2000 to January 
2013. The periodic gaps in the photometry are due to the fact 
that optical observations of Betelgeuse can only be acquired at 
Villanova from September through to early April. The optical 
photometry has good overlap with our 5 epochs of VLA/PT ra¬ 
dio data and has simultaneous coverage for 3 epochs, with the 
final two epochs only 10 days (April 2002) and thirty days (Au¬ 
gust 2003) away from the nearest photometric values. Estimates 
of the photometric values for the epochs with non-simultaneous 
coverage were obtained using the SLICK meth od within the 
PERANSO software package (I Vanmunsted|20 1 4) . The SLICK 
method subtracts the peak from the time series data and then 
carries out a Eourier transform of the residual spectrum itera¬ 
tively until a good fit is found. The red filled circles plotted in 
the bottom panel of Eigure|7]are the V magnitude values for the 


dates corresponding to our VLA/PT observations, while the sin¬ 
gle green filled circle corresponds to the date of the e-MERLIN 
observations. We simply note that there was nothing striking 
with the optical photometry around the time of the e-MERLIN 
observations and focus the remainder of this discussion on com¬ 
paring the VLA/PT data with the V band data. 

Betelgeuse is a semiregular variable with both photometric 
and radial velocity studies revealing t wo prominent peri ods; a 
short primary period of ~ 400 days (e.g.. [Smith et al.ll989h and a 
longer secondary per iod lasting ~ 2000 days (e.g.. lSanfordll933l: 
ISDencer.loneslll92^ . The long secondary period may be related 
to the convective turnover ti me of hypothetical giant convec¬ 
tion c ells on the stellar surface (IStothers & Leungll97lt[^othersl 
l2010l) . although we see no clear evidence for it in our 12 years 
of photometry. The shorter primary period has been observed 
photometrically ([p upree et al.||l987h and in radial velocity data 
([Smith et al.lll98^ iDunree et al.l T990l) and is consistently evi¬ 
dent over the 12 years of our photometry. This primary period 
displays a large range of peak-to-peak amplitude variation over 
12 years with 0.3 < AV < 0.7, which is a similar range of values 
to what were reported during the mid-19 80s, but larger tha n the 
variation reported during the early 1990s ([Krisciunasl[l992h . We 
can see from the bottom panel of Pigure|7]that our VLA/PT radio 
data (i.e., red filled circles) spans about 4 different primary peri¬ 
ods, at various phases for each of these periods. To compare the 
photometry and the radio data, we have converted the V magni¬ 
tudes to flux densities and normalized both the radio and V band 
flux densities to their values for the 2004 VLA/PT epoch. We 
note that we exclude our three 0.7 cm measurements from this 
analysis as they show no significant variability over their three 
epochs. The results are shown in the top panel of Pigure|7] It can 
be seen that there is a similar behaviour between the variability 
of the optical and the radio flux densities. 1.3 cm measurements 
are available for 5 epochs, and an increase or decrease in the ra¬ 
dio flux density is accompanied by a change of optical flux den¬ 
sity in the same direction between all pairs of epochs. The other 
radio wavelengths are measured at fewer epochs but show sim¬ 
ilar trends, although in some cases with low significance, e.g., 
between the closely-spaced second and third epochs. This shows 
that there is an irregular but distinct correlation on timescales of 
a few months. 

To understand why the optical and radio flux density might 
exhibit similar variability trends, we d i scuss the origin of the 
~ 400 day primary period. IPunree et ^ ([l987h discovered a 420 
periodic photospheric (B band; A - 453 nm) and chromospheric 
(both in the ultraviolet continuum at T = 300 nm and in the Mg ii 
h and k line emission cores) modulation of the flux density. They 
argued that this periodic modulation could not be caused by rota¬ 
tion due to the required unrealistic rapid rotation rates, nor could 
it be due to a few hot convection cells as the emergence of such 
cells would not be expected to be periodic. Periodic photospheric 
pulsations were deemed to be the most likely candidate and a 
time lag of ~ 70 days in chromospheric variation after photo¬ 
spheric variations suggested the presence of travelli ng waves in 
the atm osphere. Radial velocity observations from [Smith et al.[ 
([1989[) confirmed a primary period of 400+20 days and again as¬ 
sociated its origin to stellar pulsations. V band photometry cov¬ 
ers a large number of TiO band heads, whose formation is very 
temperature sensitive. The V band intensity change can be inter¬ 
preted to result from changes in the optical flux caused b y tem¬ 
perature changes in the photosphere ([Morgan et al.[[l99^ . How¬ 
ever, the strengths of optical TiO band heads ma y also be affected 
by tem perature effects from the chromosphere ([Lobel & Dupree! 
[2000li . Assuming the longer wavelength radio emission origi- 
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Fig. 7. Top panel: The variation of V band optical photometry with our radio flux density values for all wavelengths between 1.3 and 6.1 cm. The 
optical photometry has units of flux density and has been normalized to the October 21-30 2004 value. The radio flux densities at each wavelength 
have also been normalized to their corresponding October 21-30 2004 value. Bottom panel: The original V band optical photometry in units of V 
magnitude. The red and green filled circles are the V band photometric data at the corresponding dates of the radio observations. An example of 
the well known 420 day period is included which is believed to originate from stellar pulsations. The dashed lines in both panels should guide 
the eye. 


nates in the more extended atmosphere, one would expect a de¬ 
lay of ~ 3 - 6 years for a disturbance to propagate out from the 
photosphere to these radio formation regions , assuming a veloc¬ 
ity of 10kms“' te.g.. lO’Gorman et al.ll2012ll . Clearly, this is too 
long a timescale in comparison to the similar timescales between 
both the optical and the radio flux density variability that is evi¬ 
dent in Figure |7] 

Pulsation induced radiation field changes must be invoked 
to explain these relative quick changes. In such a scenario, pul¬ 
sation induced shocks heat the inner layers of the extended at¬ 
mosphere, resulting in a local increase in the electron density. 
Recombination (with ti mescales on the order of a few days e.g., 
iHamer & Linskvl200lh increases the global radiation field caus¬ 
ing an almost instantaneous increase in the photoionization of 
metals in the more extended atmosphere. This process can ac¬ 
count for large-scale and fast changes in the radio fluxes through 
changes in the radio opacity (i..e, kx ^ ne«ion)- We can see 
in Figure Q that our 1.3 cm flux density measurements track 
changes in the optical photometry well. We find no variation 
in the size of the extended atmosphere at this wavelength and 
so the flux density variab ility at this w a veleng th must be due to 
gas temperature changes. iHarner et al.l (12001 ) predict that most 


of the 1.3 cm emission is formed within ~ 1.8 R* and so the 
shocks might damp within this region. This is in good agree¬ 
ment with lower mass long perio d variables whose pul sation 
induced shocks damp within 2R* (iReid & MentenlflQQTh . The 
pulsation induced changes in the radio opacity at longer wave¬ 
lengths might then be the cause of the variations in the extended 
atmosphere’s effective angular radius, which is obvious in Figure 
[3 It remains unclear why the 0.7 cm flux density remains almost 
constant, but it may well be a result of the relatively large abso¬ 
lute flux density uncertainty and the possible loss of coherence at 
this short wavelength. The absence of correlated variations in the 
7 mm flux densities may also be a result of the higher densities 
and heat capacity of the emitting region, which leads to smaller 
temperature changes as compared to longer wavelengths. 


5. Conclusions 

We have used the VLA along with the Pie Town VLBA antenna 
to study the extended atmosphere of Betelgeuse at 0.7, 1.3, 2.0, 
3.5, 6.1, and 20.5 cm, over multiple epochs between 2000 and 
2004. Combining the Pie Town antenna with the most extended 
configuration of the VLA has allowed us to obtain visibilities be- 
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yond the first null at all wavelengths between 0.7 and 6.1 cm in¬ 
clusive, implying that we fully resolve the extended atmosphere 
at these wavelengths. We find that the extended atmosphere de¬ 
viates from circular symmetry at all wavelengths and that there 
may exist small pockets of gas signihcantly cooler than the mean 
global temperature of the extended atmosphere. We also find that 
the radio flux density from Betelgeuse has decreased by about 
20% in comparison to measurements from the 1970s and 1980s, 
despite the spectral index remaining almost constant at a = 1.33. 
This implies that there has been a global reduction in the gas 
temperature and/or size of the extended atmosphere. 

We hnd that both the size and temperature profile of Betel- 
geuse’s extended atmosphere as a function of wavelength is in 
good agreement with the findings of L-i-98. At 0.7 cm the ex¬ 
tended atmosphere is over twice the size of the optical photo¬ 
sphere, while at 6.1 cm it increases to about six times the size. 
We conhrm that Betelgeuse’s extended atmosphere has a low 
mean gas temperature over multiple years. The mean gas tem¬ 
perature has a typical value of 3000 K at 2R* and decreases 
to 1800K at 6R*. The temperature profile can roughly be de¬ 
scribed by the power law TgasCr) oc ® although temporal vari¬ 
ability of a few lOOK are evident at some epochs. 

The e-MERLIN 5.2 cm unresolved radio hotspots detected in 
2012 have no signature in our VLA/PT data spanning from 2000 
to 2004 even though we had sufficient spatial resolution and sen¬ 
sitivity at short wavelengths to detect them. Radio continuum 
emission is biased towards regions that are hot and ionized, the 
two most often going together, so one would expect signatures at 
multiple wavelengths. The combined constraints provided by the 
multiple epochs of spatially resolved multi-wavelength VLA/PT 
data has encouraged us to re-examine the e-MERLIN data. Erom 
this we can conclude that the e-MERLIN radio hotspots de¬ 
scribed in R-i-13 are most likely interferometric artefacts. 

We report a similar behaviour between the variability of V 
band optical photometry and the radio flux density at certain 
wavelengths. The optical variability is known to be induced 
by stellar pulsations, and so the radio flux density variability 
is likely a manifestation of stellar pulsation induced shocks. 
These shocks might damp at distances within 1.8R* and cause 
variations in the structure of Betelgeuse’s extended atmosphere 
through pulsation induced radiation held changes. The similari¬ 
ties between the variation of the optical and the radio Hux den¬ 
sities suggest that stellar pulsations could play an important role 
in exciting the extended atmospheres of red supergiants. 
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Appendix A: Deviations from a uniform disk 

In Table lA.ll we present the PT data sets where a uniform in¬ 
tensity disk plus one or more delta functions provided a better 
ht over a single uniform disk model. All of these six data sets 
contained negative point sources in the residual images after the 
best fit uniform disks from Table [1] were subtracted from their 
visibilities. 
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Table A.l. PT observations best fit by a uniform intensity disk plus one or more delta functions. 


Date 

A 

(cm) 

Shape 

Fy 

(mJy) 

^maj 

(mas) 

^min 

(mas) 

PA 

n 

Offset^ 

(mas, mas) 

Fb 

(K) 

2004 Oct 21,30 

1.3 

Disk 

16.16 + 0.37 

148+4 

106 + 4 

88 + 4 


3410+175 



Delta 

-1.21 + 0.21 




54+ 4,-1 + 5 

2585 + 355 



Delta 

-0.99 + 0.20 




-53+ 5,-3+ 7 

2735 + 345 


2.0 

Disk 

7.96 + 0.26 

173 + 8 

145+11 

126+ 18 


2490 + 235 



Delta 

-0.66 + 0.20 




48+ 13,-32 + 20 

2165 + 375 

2003 Aug 10,12 

2.0 

Disk 

6.83 + 0.34 

146+ 13 

116+ 19 

145 + 21 


3010 + 340 



Delta 

-0.85 + 0.27 




35 + 15,-25 + 24 

2395 + 330 

2002 Feb 17,18 

3.5 

Disk 

3.30 + 0.12 

225+ 11 

179 + 22 

130+ 13 


1965 + 270 



Delta 

-0.44 + 0.10 




-67+18,3 + 30 

1695 + 360 

2000 Dec 23 

0.7 

Disk 

30.00 + 0.40 

100 + 2 

90 + 2 

160 + 7 


3200+ 105 



Delta 

-0.94 + 0.25 




28 + 3,-27 + 4 

2490 + 320 


Notes, t The positional offset of the delta function from the center of the uniform disk. The first term is the offset in RA with a positive value 
denoting an offset to the east while the second term is the offset in Dec with a positive value denoting an offset to the north. 
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